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Abstract 

By using first-principles method combined with molecular dynamics simulations, the structural 
and electronic properties of mono-vacancy (MV) defect in Mobius strip formed from graphene are 
investigated. Two kinds of MV are observed depending on the local structures around defects. In 
the curved areas of Mobius strip, MV has the configuration of one pentagon and one nonagon ring 
(59-type), which is similar to that of carbon nanotube and graphene. Interestingly, the most stable 
MV appear in the twisted areas and acquires novel structure, which has two pentagon and two 
hexagon rings (5566-type) and one sp 3 hybridized carbon at the central site. Molecular dynamics 
simulations furthermore prove that in Mobius strip, both normal 59-type and novel 5566-type MV 
is a stable configuration at room temperature. Furthermore, it is expected that 5566-type MV 
is a popular MV pattern and will give rise to sp 3 -type bonding in carbon based low-dimension 
materials with chiral twisting. 
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I. INTRODUCTION 



Recently, one novel caron nanostructure, graphene based Mobius strip (GMS), has at- 
tracted extensive attention due to its special topological property, namely, with only one 
face and one edge. It is totally different from that of other carbon based low-dimension ma- 



terials, such as carbon nanotu 



synthesized in experiments 



3e and graphene. Carbon based Mobius strips were already 



2| and widely investigated by theoretical researchers. 



Their results demonstrate that GMS is a stable structure, [9[ and show novel optical 



and magnetic properties. 
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lOj Especially, it may even behave as a topological insulator. [7 1 



Since carbon based materials are usually not ideal crystals and atomic vacancies can affect 
their properties significantly, [l2- 14] such as inducing magnetism in graphene, 14] great at- 
tention has been focused on the study of atomic vacancies presented in carbon nanotube and 
graphene, 



13 



20] where atomic vacancies can be introduced unintentionally during the pro- 



cesses of synthesis or deliberately by irradiation, chemical and plasma treatments. Among 
different types of atomic vacancies in carbon nanostructures, mono-vacancy (MV) with one 
atom missing from lattice is a simple, popular and attractive one and has been identified 
clearly by experimentalists and theorists in carbon nanotube and graphene. 19M23I] In these 
systems, MV undergoes a Jahn- Teller distortion, which leads to the formation of covalent 
bond between two of three atoms located around the atomic vacancy and results in one 
five-membered ring and one nine-membered ring (59-type MV). To date, only 59- type MV 
was reported in both carbon nanotube and graphene. 

Just like carbon nanotube and graphene, to fully understand the properties of GMS for 
its potential applications, it is important and necessary to illustrate its MV properties. First 
of all, the structural configurations of the MV present in GMS should be clarified. As shown 
in our previous work,^ since each GMS is composed by two kinds of areas, curved and 
chiral twisted parts, the MV features of GMS should be investigated separately depending 
on its local structures. In this work, by taking use of first-principles method we illustrate the 
structural and electronic properties of the MV in GMS. Interestingly, besides 59-type MV 
present in the curved area, a novel 5566-type MV, which has two pentagon and two hexagon 
and acquires one sp 3 carbon atom at the central site, is observed in the twisted region. 
To the best of our knowledge, the 5566-type MV has never been reported in carbon based 
low-dimension materials before. Furthermore, due to the strain introduced by their edges, 
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graphene nanoribbons (GNR) with chiral twist were theoretically predicted. |24j-|27l| and in 




fact experimentally observed. 29[] Especially, in a very recent experiment, 30] the shapes 
of GNR with different chiral twisted degrees were clearly shown by transmission electron 
miscroscopy (TEM) images. Since graphene nanoribbons can indeed be transformed into 
carbon nanotubes through chiral twisting, [3JJ and the local configurations of the twisted 
areas of GMS are similar to that of GNR with chiral twisting, we predict that 5566-type 
MV is also observable in chirally twisted carbon nanotubes and graphene. 

II. COMPUTATIONAL DETAILS 

The calculations are performed by the widely adopted SIESTA package, in which a norm- 



conserving pseudopotential and linear combinations of atomic orbital basis sets are used. 
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The wave function is expanded with a double-^ (DZ) basis sets, and the exchange- 
correlation potential of generalized grad ient approximation (GGA) with the form of Perdew- 
Burke-Ernzerhof (PBE) is selected. [34l| The lattice vectors (50x50x50 A) are large enough 
to avoid the interactions from adjacent neighbors. All related structures are fully relaxed. 



The met 



work 



.0, 



rod and parameters used in this work are basically the same as that in our previous 
35| The accuracy of our procedure is tested by calculating the properties of MV 



in infinite graphene by taking use of 12x12 supercell. Our results are in a good agreement 



with previous ones. [23j] 



III. RESULTS AND DISCUSSIONS 



,[9 



In our previous work,|9[ the structural features and formation energies of GMS as the 
function of their width-to-length ratio are characterized. We find that one, two and three 
planar triangle regions can be observed in GMS and their formation energies increase with 
the ratio of width-to-length increasing. Typically, a GMS model with the length of 30 
armchair lines and the width of 6 zigzag lines is used to simulate the MV properties of GMS 
and such a structure has one planar triangle area. Due to its special topological structure, 
GMS can be separated to three different areas depending on local deformation. As shown 
in Fig. 1(a), A and C denote the areas with twisting, while B is the curved area similar to 
carbon nanotube. One carbon atom which is close to the center of A, B or C area is removed 
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from the GMS to create a MV and these MVs are hereafter denoted as MV-A, MV-B and 
MV-C, respectively. The formation energy (fl) of MV in GMS (labeled as Qmv (GMS)) and 
that in graphene (labeled as Qmv(G)) are respectively given by the following equations: 

Q M y(GMS) = E(GMS) - E MY (GMS) - £ atom , (1) 

^mv(G) = E{G) — E My (G) — E atom , (2) 

£atom = E(G)/N. (3) 

where E(GMS) (E(G)) and E MV (GMS) {E MV {G)) denote the total energies of 
GMS (graphene) without and with MV, respectively, and N is the number of carbon atoms 
in a 12x12 graphene supercell. The calculated results are presented in Table I, the for- 



mation energy o 



MV in graphene is 7.56 eV, which agrees well with that of the published 



results(7.4 eV[22| and 7.7 eV|23|]). The formation energies of MV-A, MV-B and MV-C are 
respectively 6.43 eV, 6.51 ev and 6.72 eV, and all of them are about 1 eV smaller than 
that of graphene (7.56 eV), indicating that compared with infinite graphene, it will be much 
easier to introduce MVs in GMS or graphene with chiral twisting due to the strain existing 
in these twisted structures. Furthermore, the most stable site for MV in GMS is located at 
its planar triangle area (A area) since the formation energy of MV-A is the smallest one, as 
shown in Table I. 

The relaxed MV structures are presented in Fig. l(b-d), and the zoomed in views are 



shown in the corresponding insets to c 
of carbon nanotube and graphene, 21 



early illustrate their structures. Similar to the MV 



23| as can be found in Fig. 1(c), MV-B shows a 
typical 59-type pattern, with one pentagon and one enneagon, and a new bond is formed 
between atom 1 and atom 2 with the bond-length of 1.59 A, which is much shorter by 0.98 
A than the distance (2.57 A) between atom and atom 1, due to the Jahn- Teller distortion. 
Interestingly, MV-A and MV-C created in the twisted areas finally evolve into a novel 5566- 
type pattern where the central carbon atom becomes sp 3 hybridized after full relaxation. 
The structural properties of MV-A and MV-C are shown in Table II, and the results show 
that the averaged distance between the central carbon and its four neighbors is about 1.64 
A, which is significantly longer than the C-C bond-length (1.42 A) in infinite graphene and 
slightly longer than that of the new bond formed in MV-B (1.59 A) and the C-C bond in 
diamond (1.55 A). Combined with the bond- angle information of MV-A and MV-C shown 
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also in Table II, we can conclude that the carbon atoms located at their central sites are of 
a sp 3 type. To further confirm this, the electronic structures are studied, and Fig. 2(a) and 
(b) respectively present the partial density of states (PDOS) of edge carbon atoms of ideal 
GMS and the carbon atoms at the central sites of ideal planar triangle area. Meanwhile, the 
PDOS of the central carbon atom of MV-A is presented in Fig. 2(c). The carbon atom at 
the edge sites of GMS is spin polarized, which is similar to the case of graphene nanoribbons 
with zigzag edge and agrees with previous work. [I], I^] Furthermore, if we compare Fig. 2 
(b) with (c), the energy gap around the Fermi level shown in Fig. 2(c) is much larger than 
that of Fig. 2(b), the PDOS of a sp 2 hybridized carbon, indicating that the central carbon 
atom of MV-A and MV-C is hybridized with its four neighbors by sj> 3 -like bonding. 

Now we investigate whehter this 5566-type MV is stable or not at room temperature. 
For this purpose, molecular dynamics simulations of MV-A and MV-B at room temperature 
(T = 300 K) are carried out in the time step of 1 fs. The corresponding results of MV-A 
and MV-B are shown in Fig. 3 (a) and (b), respectively. After running 2000 steps, the 
MV geometry of 5566-type (MV-A) and 59-type (MV-B) pattern is still kept unchanged, 
suggesting that in Mobius strip formed from graphene, not only 59-type MV but also 5566- 
type MV containing a sp 3 carbon atom is a stable configuration at room temperature. 
Finally, please note that, since the 5566-type MV can be formed at the areas with different 
degree of chiral twisting in the strip, such as A and C areas shown in Fig. 1(a), it should also 
be a popular MV configuration in carbon nanotube and graphene with chiral twists, which 
were actually observed in experiments. 29|, |30] Furthermore, we expect that the sp 3 -type 
features in this 5566-type MV can be probed by experiment, such as, X-ray photoelectron 
spectroscopy, in the chirally twisted carbon nanotube and graphene. 



IV. CONCLUSION 



The structural and electronic properties of mono-vacancy in Mobius strip formed from 
graphene are simulated. A novel mono- vacancy with 5566-type pattern with a sp 3 hybridized 
carbon atom at its central site is observed in the chiral twisted areas, while a normal 59- 
type mo no- vacancy appears in the curved area. Note that, the most stable mono- vacancy 
is the 5566 type that is located at the planar triangle area of Mobius strip. Our molecular 
dynamics simulations prove that both 5566- and 59-type mono-vacancy are stable at room 
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temperature. Furthermore, 5566-type mono-vacancy probably is a popular defects and can 
contribute sp 3 -tjpe spectral signals in the carbon based low-dimension materials with chiral 
twisting. 
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TABLE I: Formation energies of mono-vacancy. 



eV 


graphene 


MV-A 


MV-B 


MV-C 


Formation energy 


7.56 


6.43 


6.51 


6.72 
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TABLE II: The structural properties of MV-A and MV-C. Distance (d) and angle (a) are shown 
in the unite of A and degree, respectively. 





do-i 


do-2 


do-3 


d -4 


ai-0-3 


«2-0-4 




fll-0-4 


MV-A 


1.62 


1.67 


1.60 


1.68 


138.0 


147.6 


96.2 


96.1 


MV-B 


1.61 


1.69 


1.60 


1.68 


138.1 


145.7 


96.3 


96.1 
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Fig. 1 Wang.eps 




FIG. 1: (Color online) (a) The structure of ideal Mobius strip, and A, B and C are used to present 
the special areas to create MV. Red triangle shows the planar area, (b), (c) and (d) show the 
relaxed structures of MV introduced in the A(MV-A), B(MV-B) and C(MV-C) areas, respectively, 
while the enlarged structures are presented in the insets. The Arabic numbers presented in figures 
and insets are used to indicate the configurations and the atomic sites of MV, respectively. 
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Fig. 2 Wang.eps 
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FIG. 2: The partial density of states (PDOS) of carbon atom located at the edge site of ideal 
Mobius strip (a), the central site of planar triangle area of ideal Mobius strip (b) and the central 
site of MV-A (c). 
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Fig. 3 Wang.eps 
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FIG. 3: (Color online) Changes of temperature and total energy as a function of time obtained 
from molecular dynamics simulations: (a) for MV-A, and (b) for MV-B. 
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